Cooperative Ring Exchange and 
Quantum Melting of Vortex Lattices in Atomic Bose-Einstein Condensates 
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Cooperative ring-exchange is suggested as a mechanism of quantum melting of vortex lattices in a 
rapidly-rotating quasi two dimensional atomic Bose-Einstein condensate (BEC). Using an approach 
pioneered by Kivelson et al. [Phys. Rev. Lett. 56, 873 (1986)] for the fractional quantized Hall 
effect, we calculate the condition for quantum melting instability by considering large-correlated 
ring exchanges in a two-dimensional Wigner crystal of vortices in a strong 'pseudomagnetic field' 
generated by the background superfluid Bose particles. BEC may be profitably used to address 
issues of quantum melting of a pristine Wigner solid devoid of complications of real solids. 

PACS numbers: 03.75.Lm, 67.57.Fg, 73.43. Cd, 73.43.Nq 



I. INTRODUCTION 

The creation and observation of the triangular vortex 
lattices in a rapidly-rotating atomic Bose-Einstein con- 
densate (BEC) [1-3] has opened a new direction for the 
study of quantum vortex matter. Theoretical predictions 
[4] for the existence of fractional quantum Hall like states 
at even higher rotational speeds in quasi two dimensional 
atomic BEC has given a further impetus to this fascinat- 
ing field. The quantum melting of an ordered vortex 
lattice to an exotic quantum fluid of atoms at very low 
temperatures is a quantum phase transition, where one 
would like to understand the mechanism of melting and 
nature of phase transition. 

Melting of classical solids with short range inter atomic 
potential in 2D is a well studied subject, where topolog- 
ical defects play a fundamental role. In the presence of 
long range interaction, such as one component coulomb 
plasma in 2D, melting is dominated by ring exchanges [5] 
rather than topological defects. From this point of view 
the logarithmic repulsion among the imposed vortices in 
a rotating BEC provides an opportunity to study quan- 
tum melting of a 'pristine' Wigner solid with long range 
forces, that is free from the complications of solid state 
systems. 

In this letter we write down an effective Hamiltonian 
for the vortex degrees of freedom, motivated by an anal- 
ogy [6,7] between the Magnus force acting on a vortex 
moving on a two dimensional neutral superfluid fluid and 
the Lorenz force acting on a charged particle in a mag- 
netic field. We develop a theoretical approach, borrowing 
heavily from pioneering ideas of Kivelson, Kallin, Arovas 
and Schrieffer (KKAS) [8], developed in the context of 
fractional quantized Hall effect (FQHE), and suggest a 
cooperative ring exchange (CRE) mechanism for quan- 
tum melting of vortex lattices in quasi 2D atomic BEC 
and indicate a possible direction for a microscopic under- 
standing of the quantum liquid of molten vortices. 



In contrast to many recent theoretical works on atomic 
BEC which exploits an analogy between the Hamiltonian 
of a rotating neutral boson atoms and charge particle in 
an external magnetic field in two dimensions, our work 
uses the vortex (collective) coordinates directly and pro- 
vides another microscopic approach to understand quan- 
tum melting and the quantum Hall-like state that may 
be formed in these atomic system. Existing theoretical 
works focus on exact diagonalization [9] of small number 
of atoms to get some idea about quantum melting and 
the possible quantum Hall like melted states. A recent 
interesting work [10] that studies melting of vortex lat- 
tices in a rapidly- rotating 2D BEC, also shows that BEC 
is destroyed by the vortex lattices. 

Experimentally, at the present moment it is a chal- 
lenging task to produce vortex liquid state in a rapidly- 
rotating atomic BEC, in contrast to the formation of an 
incompressible liquid state of electrons in a high magnetic 
field at higher filling fractions. With the rapid advances 
in the field of laser cooled atomic gases one can antic- 
ipate to get 'snapshots' of the melted configurations of 
the vortex lattices, where CRE should leave its unique 
signatures as we mention at the end. 

In the cooperative ring exchange approach to FQHE, 
KKAS view the Laughlin quantum Hall state as a Wigner 
solid of electrons in 2D in strong magnetic field, that has 
been quantum melted by cooperative ring exchange pro- 
cesses. Briefly, ring exchange, as the name suggests, is a 
cooperative shift of a ring of contiguous particles in an 
ordered lattice (figure 1) resulting in a cyclic permuta- 
tion within the ring. While the amplitude for a quantum 
tunneling event of a specific ring of size L sites is expo- 
nentially small ~ a L (with a, the single particle tunnel- 
ing amplitude being < 1), the number of rings of size L 
is exponentially large ~ e bL . Thus the total amplitude 
~ a L e bL may exponentially diverge, if — In a < 6, lead- 
ing to a proliferation of ring exchanges and a consequent 
quantum melting. 
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This melting depends on the electronic filling fraction, 
the ratio between the density of conduction electrons and 
the density of flux quanta. For very low filling fraction, 
electrons are expected to form a Wigner crystal. At 
higher filling fraction, electrons forms an incompressible 
liquid state and exhibits quantized Hall effect. Similarly, 
we could also expect the quantum melting of the vortex 
lattices depends on the vortex filling factor, the ratio of 
the total number of vortex to the total number of boson. 



II. HAMILTONIAN OF THE VORTICES IN A 
ROTATING QUASI-2D BEC 

We consider a large number of vortices in a rapidly ro- 
tating quasi-2Z? BEC; the condition for the condensate in 
a quasi-2D trap is p = p n g 2 < Tilo z and the atoms in the 
condensate should be in the atomic lowest Landau level 
generated by the fast external rotation is p g 2 < 2htt. 
Here, po is the boson density, g 2 = 2^/2~nfiuj z a z a is the ef- 
fective interaction strength in quasi-2D Bosc system [11], 

where a z = 



N v 



5— is harmonic oscillator length along 

the z-direction with oj z is trap frequency in the axial (z) 
direction and £1 is the trap rotational frequency. Also, 
a is the s-wave atomic scattering length. A vortex in a 
fluid is an excitation in which each fluid particle is given 
an angular momentum to relative to the vortex center. 
Here, we treat a vortex as a point particle moving under 
the influence of the Magnus force. The Magnus force is 
an effective interaction between superfluid particles and 
vortices in relative motion [6,7]. The force acting on a 
single vortex [6] is then 



F = v x z(2ttT%pq). 



Here, v is the vortex velocity relative to the superfluid 
particles and po is the superfluid particle density. The 
Magnus force is equivalent to the Lorentz force acting 
on a charge particle (e) in a magnetic field. Hence, 
eB e ff = 2-K%p is the pseudo magnetic field. 

The interaction potential between two vortices sepa- 
rated by a distance r is 



2^po/r 
m U 



where £ 



' is the coherence length of the vortex 

nan ° 



core and m is the mass of a superfluid particle. The 
above potential is valid only when the distances between 
two vortices is greater than the coherence length. No- 
tice that the interaction strength between two vortices 
depends on the superfluid density as well as the s-wave 
scattering length a. 

The Hamiltonian of a rotating BEC containing vortices 
can be written in terms of center of vortices (collective 
coordinate) as [6] 



(Pi - irhpoz x r l ) 2 2i:h 2 p 



2m v 



mo 



i<j 



In 



(3) 



where N v is the total number of vortices. The effec- 
tive vortex mass m v = Trpo£ 2 can be in principle derived 
from a microscopic approach [6]. Since the coherence 
length is very small, the vortex mass also becomes small. 
This Hamiltonian is similar to that of a charged parti- 
cles moving under the influence of the Lorentz force by 
a magnetic field B c g. The pseudo vector potential due 



r x B c ff. For N v num- 



to the Magnus force is A c gf = 
ber of vortices in an area A, one gets the vortex filling 
factor, v v = ^j-jjir^ = ^r- Notice that the vortex fill- 
ing factor v v is just inverse of the bosonic filing factor 
Vb = j!-. For large N v the vortex density is approxi- 
mately uniform and N v = 2m £ A . N is the number of 
the superfluid particles. The effective magnetic length 
is l " = \J = ■ Thc pscudomagnetic field gen- 

erated by the background superfluid particles leads to 
the quantization of the cyclotron motion and producing 
Landau levels of the vortices. The eigen spectrum of 
the single vortex Hamiltonian are uniformly spaced with 
energy gap Tilu c h, where oj c s = 2 ^ Po is the effective cy- 
clotron frequency. The limit of m v — ► and/or large 
superfluid density (po) is equivalent to the vortices are 
in the lowest Landau level (LLL). We can project thc 
Hamiltonian onto the LLL and the corresponding nor- 
malized wave functions are degenerate eigenfunctions of 
the angular momentum m is, 



(!) V(z) 



m+lr 



TO 



0,1,2,.... (4) 



where z = x + iy and (x, y) are the position coordinate 
of a vortex. When the vortices are confined to the lowest 
Landau level, (i.e. the cyclotron degrees of freedom are 
confined to the LLL), the kinetic degrees of freedom of 
the vortices are frozen, since the spacing between Landau 
levels, Jiuieff, is large compared with all other energies in 
the problem. The remaining degrees of freedom are the 



(2) vortex guiding-center coordinates, R 



fc(P x 



The guiding center coordinate R specify the center of 
a Gaussian- localized probability amplitude of width Iq. 
These coordinates have no kinetic energy. Hence the 
vortices in thc LLL will remain localized about a given 
guiding center coordinate R indefinitely. 



III. COHERENT STATE PATH INTEGRATION 

In this section, we would like to review the coherent 
state path integral formalism. For detailed derivations, 
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please consult the references [12,13]. In symmetric gauge, 
the wave function of a vortex in the LLL with guiding- 
center position R is 



Mr) 



exp 



R| : 



All 



+ 



i(r x R).z 



2l 2 



(5) 



It has the same form as a coherent state in a two- 
dimensional phase space [12]. Here, the state label R 
is a continuous variable. The coherent state overlap is 
given by, 



< R1IR2 >= exp 



|Ri — R2 



+ 



i(Ri x R 2 ).z) 



2Zg 



(6) 



This coherent state |R > forms a nonorthogonal, over- 
complete basis. Nevertheless, the projection operator P 
onto the LLL is given by, 



d 2 R 
2irl 2 



R >< R| 



(7) 



which is unity within the LLL since < Ri|P|R2 >=< 
Ri|R 2 > 

We use the coherent state path integral [12,13] expres- 
sion for the partition function to calculate the tunneling 
coefficient of a vortex. The partition function for 2D in- 
teracting vortices in a pseudomagnetic field due to the 
Magnus force is 



Z{v v )=Tre- pHv . 



(8) 



Here, we discuss the main features of this formalism for a 
single vortex in the LLL in the complex plane. This can 
be generalized for many vortex system very easily. The 
coherent state in the complex plane is 



\R >= 



1 



exp 



zR* - z*R 



All 



exp 



R\ 
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, (9) 



where R = X + iY is the guiding center coordinate of a 
vortex in the complex plane and the asterisk denotes the 
complex conjugation. The coherent state overlap in the 
complex plane is 



< Rj\Rk > — exp 



exp 



{ Al 2 ' 



\Rj - R k \ 



2/ 2 



<««,>=/n^n<* 



. j+1 \e eV \Rj >, 



(12) 



where i?o = Ri,R n +i — Rf. In general, the Hamiltonian 
can be written 



V(R*,R k ) 



< R 3 \V\R k > 
< RARk > 



(13) 



The matrix element can be written as, 



<Rf\e-> H \Ri> = /n^n<^ + ^> 

x [l- l eV(R* +1 ,R 1 ) + 0(e 2 )] (14) 

We are neglecting terms of 0{e 2 ) and higher order terms 
by standard procedure. Using Eq.(lO), and 



dRj 



R j+1 - Rj 



we obtain 



Z(y v ) 



nd 2 Rk 
7^J^ ex P 



' 51 47? S 3 



1=0 



where 

s > = { R *i^r ~ d ~it Rj+1 ) V{R ^ R 

The above path integral can be written as 



(15) 



(16) 



J D[TL]t 



iS[R] 



where 



Jo 



dt 



4^Kf ~ d J R ^^' 



(17) 



(18) 



(19) 



[m{Rj - R k ) - (R* - R* k )Rk 



Al 2 



(10) 
(11) 



Now the path integral representation of the partition 
function Z(v v ) can be obtained in the usual way. First, 
we split the inverse temperature (5 into a large number of 



equal intervals e = (3/n, i.e., e~^ v is written as [e _ey ] n , 
and then insert the projection operator P at each in- 
finitesimally small interval. Then, 



This action is linear in time derivatives and hence dis- 
continuous paths have finite action. It implies that the 
coherent state path integral is dominated by discontin- 
uous paths and the limits is ill defined. Despite these 
difficulties, the continuum version of the path integral 
can be used to develop a saddle-point approximation for 
the partition function [8] . We are interested in the semi- 
classical limit when V(R) is a slowly varying function of 
its argument over the length scale lo and we can use the 
saddle point approximation to evaluate the path integral. 

The single vortex path integral can be generalized to 
many vortex path integral directly. The action for many 
vortex is 
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S(K)= f dt ^(R J xR 3 )i + ^y(R,-R,) 

J ° [ 3=1 j<k 

(20) 



where 



V(R) = (Mr)\V(r)\<fa.(r)) 



(21) 



is two-body interaction potential in coherent states rep- 
resentation. 

In the saddle point-approximation, the classical path 



is obtained by minimizing the action, SR ^ |fl=_R e = 0. 
The classical paths satisfy the following equations of mo- 
tion 



x z, where Vj 



J2y 

k^3 



;R,-R fe ). (22) 



The path integral can be expressed as a sum over sad- 
dle point contribution in which the contribution of paths 
in the neighborhood of each classical path is evaluated by 
expanding the action to quadratic order in R — R c . The 
partition function Z is calculated within the semiclassical 
approximation. The partition function can be expressed 
as a sum over classical paths, assuming the vortices to be 
bosons, Z = Y,c D [ R c] e ~ S[Rc \ wherc D [ R c\ is the fluc- 
tuation determinant and S[R C ] is the action evaluated 
along the classical path. There are interesting issues [14] 
about the statistics of vortices in a compressible super- 
fluid such as ours, and in the present paper we assume 
them to be bosons. The partition function Z can be writ- 
ten by considering only the leading order contribution as 
[8] Z = Z G Y, c D[Rc]e- Sa[R "\ where S = S + S. Let us 
consider the contribution of a single large exchange ring 
to Z. The real part of the action would be a$L (see the 
sec. V), where L is the number of vortices in the ring 
and «o is independent of path. The fluctuation deter- 
minant [8] is D[R C ] ~ exp[— 5aL + O (In L)]. Here, 5a 
is a real constant which renormalize cto- The imaginary 
part, the phase change as a cooperative motion along a 
ring is 8 = § eA ff.cJl = 2irN, where N is the number of 
the supcrfluid particles enclosed by the ring [14]. This 
is the analog of the Bohm-Aharanov phase factor for a 
charged particle moving in a magnetic field. So the par- 
tition function becomes, Z ~ exp[— ah ± i2irN], where 
a = ao + 5a. 



IV. COOPERATIVE RING EXCHANGE 
MECHANISM 

How does the vortex lattice melt? To understand the 
melting of the vortex lattices, the Lindemann criteria can 
not be used here since it is used in the melting of clas- 
sical solids. The vortices are not executing almost in- 
dependent thermal motions as in a classical solid. The 



dynamics of the present problem is governed by a Hamil- 
tonian with only first-order time derivatives, which give 
rise to its own peculiar properties. If we consider a rigid 
Wigncr solid and allow one ring of vortices to tunnel 
coherently they see a periodic potential with the period- 
icity of the lattice (Fig 1). If we observe the coherent 
motion of one chain over a long time compared to the 
tunneling time (to), the potential that it sees will not 
be periodic. The physically important rings being one- 
dimensional and long, this can result in the destruction of 
the long-range order along the chain rather easily. This 
in turn will feed back and affect the rest of the neigh- 
borhood, resulting possibly in a molten state. This will 
also result in the path of the wave packets of vortices 
being displaced away from the edges of the triangle of 
the lattice. This means that the self-consistent potential 
seen by a vortex no longer has a component which has 
long-range order. 
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FIG. 1. A schematic diagram of cooperative ring exchange 
events on a ring and line. The vortices are indicated by the 
small circles. 



V. CALCULATION OF THE TUNNELING 
COEFFICIENT 



The numerical value of the tunneling coefficient a(v v ) 
determines whether the vortices form the liquid state or 
Wigner crystal. To estimate this tunneling coefficient 
we consider the following simple exchange path which 
is shown in figure 1. Consider the path in which one 
row of vortices exchanges one step in the X direction 
in the background of the static potential of all other 
vortic es, Xj (0) = X t (0) + d and Y t {J3) = F;(0), where 

d = y/ lo is the lattice constant of the Wigner crys- 
tal for a given density v v . There is no net phase changes 
since this straight path does not enclose any area. We 
are imposing the periodic boundary conditions in the X 
direction, X 1 (t) = X i+L (r). 

For \Yj \ <§C d, the two-body interaction potential in the 
coherent state representation (given in the eq. 21) can 
be approximated by 



4 



v_ 

m 



^{E[fn ! + ^d-cos(^))j (23) 
+ I Y^ K xi3 - k )( X 3 - X kf + K y {] - k){Yj - Y k ) 2 ]}, P laced b y (j-Wxt, but £\ i 2 ^(i) is diverging linearly 



evaluate the fluctuation determinant we have to take the 
continuum limit of the effective action S c s. To take a 
continuum limit of the effective action, (<f>j — (f>k) is rc- 



j>k 

where X'jS and Y's are in units of the lattice constant d. 



Here, K x (j-k) = a d x Jix k W and K y (j-k) = £0t\r c 
are evaluated along the classical paths. The best fit to 
the actual potential is obtained with Q x /Q y ~ 0.6 for 
v v ~ 1/2 and Q x , Q y and K x (j), K y (j) are weakly depen- 
dent on v v . The calculation of the fitting parameters Q x 
and Q y is given in the Appendix. Notice that Q x /Q y < 1 
implies that when one-dimensional chain moves coher- 
ently in the X direction, the potential barrier is much 
less than in the X direction compared to that of the Y 
direction. 

The dimensionless Euclidean action is given by 

s/ZvJo L Y 2 3 3 2 1 
+ ^(l-cos(2^))] 

+ \Y,\. K *V ~ k ^ X i - + K vU - k )( Y i - Y k)% 

j>k 

Since S is a quadratic form in Yj, the motion in the Y 
direction can be integrated out exactly. After doing the 
Yj integration, we get an effective action S e s for the X 
motion, with a quadratic kinetic energy, 

S cff = - J— f dr[J2[UjM(j ~ k)4> k (24) 
V3irv v Jo ^2 

+ \[K x {j - k)(<f>j - 4> k f] + E T^ 1 - cos(^))], 

3 

where r is the imaginary time variable, (M(j — k))^ 1 = 

|[(% + E, K y (j)Sj k )-E j<k K v (j-k)} and ^ = 2ttX 3 . 
S c s is the effective action for a one-dimensional sine- 
Gordan chain. The classical path satisfying the boundary 
conditions are <f>j(0) — and <f>j((3) — 2tt, corresponds to 
the simultaneous coherent motion of all the vortices, i.e. 
4>j{r) = 4>o( T ) + 27rj. Due to the simultaneous coher- 
ent motion of all the vortices, the above effective action 
becomes, 



d 2 v„ 



1 



'off 



+ Qx{l ~ cos(^)) 



• (25) 



By using the Euler-Lagrange equation of motion, one 
can calculate the <j>. Hence, the above effective action 
along the classical path becomes, S[R C ] = a L, where 

doifv) — ~tt — \/rf L - The a is independent of K's. To 

V o7rv v y Wy 



since K x (j) ~ 1/j 2 . This is an infrared divergence and 
the continuum model must be constructed by taking the 
upper cut-off limit carefully. Here, we do not calculate 
the 5a which is a non-trivial task. 

Kivclson et al. [8] has given an extensive discussion of 
how to map Wigner crystal of electrons in a magnetic 
field into the discrete Gaussian model. Following the ref. 
[8], one can map the sum over all classical paths to a 
sun over classical spin configurations. All the contribu- 
tions of ring exchanges happening in a time interval To 
are summed by modeling the change in the action by a 
discrete Gaussian model in an imaginary field [15], 



H 



DG 



Oto{Vv) 



E 

<A : 7> 



{Sx-S^f+ihiv^Sx, (26) 



where < A, 7 > denotes a nearest-neighbor pair on the 
dual lattice and Sa is an integer variable associated to 
every triangle in the lattice. «Sa counts the number of 
clockwise minus counterclockwise ring exchanges that 
surround a plaquette A. The function ao{v v ) is a mea- 
sure of the tunneling barrier. The function h{y v ) is the 
phase factor which arises as a result of the pseudo mag- 
netic flux enclosed by the exchange rings. This model is 
known to have a phase transition [15] at a critical value 
of a = a c {v v ) <~ 1.1 [8]. For a{v v ) > a c {v v ), the ground 
state is a vortex Wigner crystal and for a(v v ) < a c (v v ), 
the ground state is a quantum mechanical vortex liquid 
state. In our calculation we find that the quantum melt- 
ing will occur at v v ~ i. The current experiments [1,2] 
with f, C ^ are in the regimes of vortex lattice ground 
state. So our result is consistent with the experimental 
results, but it does not match very well with the other 
theoretical results [9,10]. The numerical calculation is 
based on the small number of atomic bosons as well as 
small number of vortices. In our approach we assumed a 
large number of atomic bosons and vortices. We believe 
that this discrepancies is related to the system size. One 
would say that a large system has been considered in the 
ref. [10] and calculated the melting condition which is 
comparable to the numerical result [9]. In the ref. [10], 
first they have calculated the root mean square of the 
displacement from the equilibrium position of a vortex 
in terms of the filling factor v v . Then they have used 
the Lindemenn criterion and assumed that the melting 
will occur when the fluctuation of the vortex position is 
0.15d to get the melting condition which is close to the 
numerical result [9]. Although the Lindemann criterion 
gives a reasonable description of the melting of a classical 
solid, there is little evidence that it can be applied to the 
melting of a vortex lattice. The vortices are intrinsically 
quantum objects whose equation of motion are quite dif- 
ferent from those of atoms in a harmonic crystal. 
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VI. SUMMARY 

In this paper, we treated the vortices as a new de- 
grees of freedom and considered a model Hamiltonian 
of interacting vortices. Later, we assumed the vortices 
arc in the lowest Landau level due to the low mass of 
the vortices and the high densities of the supcrfiuid Bosc 
particles. The concept of cooperative ring exchange is 
introduced to explain the mechanism of quantum melt- 
ing of the Wigner crystals. Finally, we estimated the 
tunneling coefficient which determines the condition for 
quantum melting instability of the vortex lattices. Lat- 
est experiments [1] with N ~ 10 5 , N v ~ 10 (y v ~ 1(T 4 ) 
and [2] with N ~ 10 7 , N v ~ 100 {u v ~ 10~ 5 ) are in the 
regime in which the ground state is a vortex lattice. It 
is a challenge for experimentalists to produce a vortex 
liquid state in a rotating Bose condensed state. 

Our present work, resulting in a discrete Gaussian 
model (equation 10) predicts Laughlin like even denomi- 
nator bosonic vortex filling fraction v v = |, to emerge on 
quantum melting. We can also determine the asymptotic 
form of the wave functions [16]. Along with a rich phase 
structure the discrete Gaussian model also determines 
the nature of the quantum melting transition. To the ex- 
tent the vortex degrees of freedom retain their identity, 
the results of CRE approach may remain valid in the 
quantum melted region. This needs to be investigated 
further. 

As mentioned earlier CRE processes should leave its 
finger print as specific fluctuation patterns (figure 1) that 
preempts quantum melting. It should be interesting to 
look for snapshots of such displaced large rings in the 
actual vortex lattice imaging. 

APPENDIX A: CALCULATION OF THE 
PARAMETERS Q x AND Q Y 

Here, we describe how to calculate the parameters Q x 
and Q y . We consider the simplest possible exchange 
path, namely one line of vortices shifting coherently 
within the Wigner crystal. When the line C is displaced, 
we have = T, + &b~iec, with Si^c unity if and only 
if lattice site i lies on the line in question. The matrix 
element of the potential between two vortices in coherent 
basis state is 

V(K) = (d> R (r)\V(r)\Mr).) (Al) 

Accordingly, the energy of the displaced line configura- 
tion relative to that of the perfect Wigner crystal is 

This sum can be broken up into three terms. The first 
term includes all pairs (i, j) in which both sites i and j lie 



off the line. This contribution to AE is zero. The second 
term involves all pairs (i, j) where one of the sites, say i, 
is on the line and the other, j, is off the line: 

AE 2 = ]T [y(T,+d-T,)-V(T 4 -T,)]. (A3) 

Clearly the line energy is extensive, hence the energy per 
tunneling of the vortex can be written 

U(d) = AE 2 /L = J2 lV(Tj - a) - V(Tj)] , (A4) 

where we have chosen the origin to lie on the line. The 
third and final term is that arising from both i and j on 
the line. Since the tunneling is cooperative, this contri- 
bution to the classical action vanishes. 

By allowing one line of vortices to tunnel coherently 
along the line, one can fit the change in energy AE into 
a periodic potential with the appropriate choice of the 
parameter Q x . On the other hand, by allowing one line 
of vortices to tunnel coherently perpendicular to the line, 
one can fit the change in energy into a quadratic potential 
with appropriate choice of the parameter Q y . 
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